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SUMM  AR  Y(U) 

A  recent  semi-empirical  model  of  near-surface  bubble  size  distributions  is 
compared  with  the  measurements  of  Walsh  and  Mulhearn  (1987),  and  is  shown  to 
provide  a  satisfactory  description,  for  bubble  radii  greater  than  approximately 
70  tun.  A  strength  of  the  method  is  that  it  can  relate  whitecap  coverage,  which  is 
readily  observable,  to  bubble  numbers,  and  the  former  is  very  likely  to  be  a  better 
parameter  than  wind-speed,  the  normally  used  variable,  for  estimating  bubble 
numbers. 
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1.  INTRODUCTION 

Bubbles  from  breaking  surface  waves  have  a  number  of  effects  on  acoustic  propagation.  They  absorb  and 
scatter  sound  and  also  refract  it  by  altering  the  sound-speed  profile.  Absorption  and  scattering  can  become 
significant  at  frequencies  above  approximately  20  kHz.  Bubbles  also  have  a  role  in  near-surface  ambient 
noise,  and  an  influence  on  the  onset  of  hydrodynamic  and  hydoacousfic  cavitation.  In  addition  they  have 
a  strong  influence  on  particulate  matter  transfer  to  the  sea-surface,  on  gas  transfer  across  the  sea-surface, 
and  on  the  generation  of  marine  aerosols. 

In  this  report  the  simple,  semi-empirical  model  of  Monahan  (1988)  for  bubble  size  distributions  just  below 
the  sea  surface  is  compared  with  the  experimental  results  of  Walsh  and  Mulhearn  (1987)  and  is  found  to 
provide  a  satisfactory  explanation.  There  is  no  depth  dependence  in  the  model.  Walsh  and 
Mulhearn  (1987)  obtained  bubble  size  distributions  with  in-sifu  photographic  measurements,  at  sea. 


2.  BUBBLE  SIZE  DISTRIBUTIONS 

Monahan  (1988)  presented  a  model  for  near-surface  bubble  size  spectra  which  related  the  size  spectrum  to 
white<ap  coverage,  W,  on  the  basis  of  a  simple  bubble-plume  model  which  "worked  backwards"  from  an 
observed  relationship  for  the  rate  of  marine  aerosol  production.  This  model  is  valid  for  "bubbles  that 
reach  the  sea  surface,  burst,  and  produce  spray  droplets"  (quoting  Monahan  (1988)).  Monahan  showed 
good  agreement  with  experiments  of  /ohnson  and  Cooke  (1979),  Kolovayev  (1976),  and  Baldy  and 
Bourgel  (1985).  According  to  the  model,  if  N(r,o)  is  the  near-surface  concentration  of  bubbles  with  radius 
r,  just  below  the  sea-surface,  then 


N(rp)  =  v(r)-’rfP(r)H(r)W(U) 


(1) 


where  v 
/ 

P 

H 

W 

U 


a  bubble's  terminal  rise  velocity, 

number  of  jet-drops  produced  per  bursting  bubble  at  the  sea-surface, 

(from  observations  1  <  j  <  5), 

the  jet-drop  fraction  of  the  bubble-generated  aerosol, 

bubble  flux  to  the  surface  within  an  active  white-cap, 

irtstantaneous  fraction  of  the  sea-surface  covered  by  active  white-caps,  and 

wind  speed  at  10  m  elevation  (ms*’). 


V  is  a  function  of  bubble  radius,  but  also  depends  on  whether  or  not  the  bubbles  are  hydrodynamically 
clean  or  dirty.  From  Thorpe  (1982) 


V  =  2/9(r^g/t>)[(y7 +2y)'/^  -  y  1,  for  dirtj  bubbles. 


where  g 
u 

y 


acceleration  due  to  gravity, 
kinematic  viscosity,  and 
10.82  tP/grJ. 
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While  for  clean  bubbles 


V 


1  /3  g/u,  for  Reynolds  Number 


2r  v/u  <  0.33, 


and 


V  =  r^g/dl- 2/11 +(1+ 0.091  gi3/u)l)18u) 


tor  larger  Re  such  that  (Re)  '''*  <<  1. 
P,  H  and  W  are  given  by  Monahan  as: 


P(r)  =  1- 0.715  e-0-030r. 


(2) 


H(r)  =  3.13xlO«{(l +0.057A>  05)xlol  '’exp(-BxB)}  /a 


(3) 


where  A  =  (8.77  x  lO'^  r  +  0.98), 

B  =  (0,380 -log  (8.77  xl0-2r  +  0.98))/0.65, 

and 


W  =  3,84  xl(r‘U34i. 


(4) 


Monahan  obtained  the  expressions  for  P,  H  and  W  from  Woolf  et  al  (1987),  Monahan  (1965)  and  Monahan 
and  O'Muircheartaigh  (1980),  respectively.  The  predicted  range  of  bubble  concentrations  for  W  =  1  (i.e.  a 
completely  whitecap  covered  sea)  is  shown  by  the  shaded  region  of  figure  1,  where  the  upper  bound  w:is 
calculated  using  Thorpe  (1982)'s  rise  velocity  for  dirty  bubbles  and  )  =  1,  and  the  lower  bound  :xas 
calculated  using  the  rise  velocity  for  clean  bubbles  and  J  =  5.  The  slope  of  Monahan’s  curves  indie  .tes  an 
r-4.3  dependance,  approximately.  During  a  number  of  the  experimental  runs  of  Walsh  and  Mulheam 
(1987)  to  obtain  in-situ  photographic  measurements  of  bubble  size  spectra  at  sea,  some  photographs  were 
obtained  which  showed  an  extremely  large  number  of  bubbles.  These  were  assumed  to  have  been  taken 
very  close  to  a  breaking  wave,  within  a  bubble  "cloud".  The  experimental  points  on  figur?  1  are  bubble  size 
spectra  from  some  of  these  "clouds".  It  can  be  seen  that  for  two  of  them  agreement  w  th  Monahan's  model 
for  W  =  1,  is  satisfactory  at  least  for  larger  radii.  It  can  also  be  seen  that  the  ir  xlel  provides  an  upper 
bound,  as  it  should,  for  W  ==  1,  if  the  model  is  valid.  v 
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Bubble-size  spectra  (averaged  over  the  32  photographs  obtained  in  any  one  experimental  run  of  17  min 
duration)  are  presented  in  figures  2,  3, 4  and  5  for  wind  speeds  of  14  ms"' 11  to  13  ms"',  8  to  10  ms"' and  6  to 
7  ms"’  respectively.  (From  equation  (4)  wind  speeds  of  14  ms"’,  11  ms"’,  8  ms"’  and  6  ms"’  correspond  to 
W's  of  3%,  1.4%,  0.5%  and  0.2%,  respectively.)  It  can  be  seen  that  the  experimental  points  scatter  about 
Monahan's  model  curves.  On  each  of  figures  3  to  5  the  experimental  distribution  averaged  over  the 
appropriate  wind-speed  interval  is  presented  and  reasonable  agreement  with  Monahan's  model  is 
obtained.  (One  recent  experiment  by  Tate  (1987)  appears  to  disagree  with  these  findings  in  that  bubble 
numbers  considerably  exceed  those  predicted,  however  the  sea  state  encountered  was  much  rougher  than 
expected  given  the  measured  wind  speed.) 


3.  DISCUSSION 

The  main  weakness  in  the  model  is  equation  (4),  which  relates  white-cap  coverage,  W,  solely  to 
wind-speed.  W  does  depend  on  other  factors  and  it  is  possible  that  if  direct  measurements  of  W  had  been 
obtained  in  the  experiments  of  Walsh  and  Mulheam  (1987)  better  agreement  would  have  been  reached 
and  less  scatter  obtained  in  figures  2  to  5.  The  regressions  of  bubble  numbers  against  wind-speed  (Walsh 
and  Mulhearn  (1987)  found  numbers  U^  ^)  and  of  W  against  wind-speed  (see  equation  (4))  show 
considerable  scatter  but  one  would  expect  white-cap  coverage  and  bubble  numbers  to  be  closely  related  so 
that  W  would  be  a  superior  parameter  for  estimating  bubble  populations.  It  is  also  more  directly 
obtainable  than  wind-speed  by  remote  sensing  techniques. 

Figures  1  to  5  also  demonstrate  that  the  numbers  of  bubbles  with  radii  less  than  approximately  70  mm 
are  less  than  Monahan's  model  predicts.  The  model  is  only  suitable  for  bubbles  which  reach  the 
sea-surface,  burst  and  produce  jet-drops.  Smaller  bubbles,  which  are  lost  by  dissolution  or  by  coalescence, 
before  they  reach  the  surface,  are  not  accounted  for.  There  are  also  unexplained  differences  between  the 
radii  of  the  maxima  in  the  bubble  spectra  of  Walsh  and  Mulheam  (1987),  Johnson  and  Cooke  (1979)  and 
Kolovayev  (1976).  Flowever,  Monahan's  model  satisfactorily  describes  the  bubble  size  spectra  at  larger 
radii. 

Monahan  (1988)  does  not  treat  the  depth  dependence  of  bubble  size  distributions.  Due  to  the  smallness  of 
the  current  data  base  it  is  not  possible  to  mathematically  specify  this  too  exactly.  The  latest  attempts 
are  those  of  Wu  (1988)  and  Hall  (1989).  The  strength  of  Monahan's  model  is  that  it  can  relate  bubble  size 
distributions  to  a  simple,  readily  observable  parameter  -  whitecap  coverage.  This  is  more  satisfactory 
than  using  wind  speed  given  the  complex  relationship  between  wind  speed  and  sea-state.  Once  the  near 
surface  distribution  is  specified  empirical  curve  fits  or  turbulent  dispersion  models  can  then  be  used  to  find 
the  distribution  of  bubbles  with  depth. 


4.  CONCLUSION 

For  bubble  radii  larger  than  70  pm  Monahan's  model  provides  a  satisfactory  description  of  near-surface 
bubble  size  spectra.  Without  coincident  measurements  of  bubbles  and  whitecap  coverage  it  is  not  possible 
to  definitely  say  that  the  latter  provides  a  better  indicator  of  bubble  numbers  than  wind-speed,  but  this 
does  seem  very  likely. 
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Figure  1.  Bubble  Cloud  size  distributions  compared  with  Monahan's  model  for  W  =  1 


/ 


WSRL-TM-55/88 


\\SRL-TM-55/88 


8 


Figure  3.  Bubble  size  distributions  for  11  to  13  ms"’  compared  with  Monahan's  model 
for  12  ms"', _ ,  average  of  the  3  measured  distributions 
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Figure  5.  Bubble  size  distributions  for  6  to  7  tns''  compared  with  Monahan's  model 
for  6.5  ms"', . ,  average  of  the  5  measured  distributions 
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